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Abstract. A parameter estimation problem is considered for a diagonaliazable 
stochastic evolution equation using a finite number of the Fourier coefficients of 
the solution. The equation is driven by additive noise that is white in space and 
fractional in time with the Hurst parameter H > 1/2. The objective is to study 
asymptotic properties of the maximum likelihood estimator as the number of the 
Fourier coefficients increases. A necessary and sufficient condition for consistency 
and asymptotic normality is presented in terms of the eigenvalues of the operators 
in the equation. 

1. Introduction 

In the classical statistical estimation problem, the starting point is a family P e of 
>- ■ probability measures depending on the parameter 9 in some subset O of a finite- 

dimensional Euclidean space. Each P e is the distribution of a random element. It is 
assumed that a realization of one random element corresponding to one value 6 = 6q 
of the parameter is observed, and the objective is to estimate the values of this 
parameter from the observations. 

o; 

The intuition is to select the value 6 corresponding to the random element that is 
most likely to produce the observations. A rigorous mathematical implementation of 
this idea leads to the notion of the regular statistical model [3] : the statistical model 
(or estimation problem) P e , 9 6 9, is called regular, if the following two conditions 
are satisfied: 

• there exists a probability measure Q such that all measures P e are absolutely 
continuous with respect to Q; 

• the density dP /dQ, called the likelihood ratio, has a special property, called 
local asymptotic normality. 

If at least one of the above conditions is violated, the problem is called singular. 
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In regular models, the estimator 9 of the unknown parameter is constructed by max- 
imizing the likelihood ratio and is called the maximum likelihood estimator (MLE). 
Since, as a rule, 9 ^ 9o, the consistency of the estimator is studied, that is, the con- 
vergence of 9 to 9q as more and more information becomes available. In all known 
regular statistical problems, the amount of information can be increased in one of two 
ways: (a) increasing the sample size, for example, the observation time interval (large 
sample asymptotic); (b) reducing the amplitude of noise (small noise asymptotic). 

In finite-dimensional models, the only way to increase the sample size is to increase 
the observation time. In infinite-dimensional models, in particular, those provided by 
stochastic partial differential equations (SPDEs), another possibility is to increase the 
dimension of the spatial projection of the observations. Thus, a consistent estimator 
can be possible on a finite time interval with fixed noise intensity. This possibility 
was first suggested by Huebner at al. [2] for parabolic equations driven by additive 
space-time white noise, and was further investigated by Huebner and Rozovskii [3], 
where a necessary and sufficient condition for the existence of a consistent estimator 
was stated in terms of the orders of the operators in the equation. 

The objective of the current paper is to extend the model from [3] to parabolic 
equations in which the time component of the noise is fractional with the Hurst 
parameter H > 1/2. More specifically, we consider an abstract evolution equation 



u{t) + f (A + 9Ai)u{s)ds = W H (t), (1.1) 

Jo 



where Aq, A\ are known linear operators and 9 G C M. is the unknown parameter; 
the zero initial condition is taken to simplify the presentation. The noise W H (t) is a 
formal series 

oo 

w H (t) = j2 w f® h J> ( L2 ) 

where {w? , j > 1} are independent fractional Brownian motions with the same Hurst 
parameter H > 1/2 and {hj, j > 1} is an orthonormal basis in a Hilbert space H; 
H = 1/2 corresponds to the usual space-time white noise. Existence and uniqueness 
of the solution for such equations are well-known for all H G (0, 1) (see, for example, 
Tindel et al. [T3J Theorem 1]). 

The main additional assumption about (11.11) . both in [3] and in the current paper, is 
that the equation is diagonalizable: {hj, j > 1} from (11.21) is a common system of 
eigenfunction of the operators Aq and A\. 

Aohj = Pjhj, Ai = Vjhj. (1.3) 

Under certain conditions on the numbers pj, Vj, the solution of (11.11) is a convergent 
Fourier series u(t) = X^>i u j{t)hj, and each Uj(t) is a fractional Ornstein-Uhlenbeck 
(OU) process. An iV-dimensional projection of the solution is then an iV-dimensional 
fractional OU process with independent components. A Girsanov-type formula (for 
example, from Kleptsyna et al. [7, Theorem 3]) leads to a maximum likelihood 
estimator 9n of 9 based on the first N Fourier coefficients uu . . . , un of the solution 
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of ( 11. ip . An explicit expression for this estimator exists but requires a number of 
additional notations; see formula (13.81) on page [H] below. 

The following is the main results of the paper. 

Theorem 1.1. Define fij = Qvj + pj and assume that the series £) 7 -(l + \(J>j\) ~ 7 

converges for some 7 > 0. Then the maximum likelihood estimator On of is strongly 
consistent and asymptotically normal, as N — > 00, if and only if the series ^\. rfnj 1 
diverges; the rate of convergence of the estimator is given by the square root of the 

( \ 1 / 2 ~ 

partial sums of this series: as N — > 00, the sequence (^2j < N u jt l J 1 ) (On — #) 

converges in distribution to a standard Gaussian random variable. 

If the operators Aq and A\ are elliptic of orders m$ and mi on L 2 (M), where M is a 
(i-dimensional manifold, and 1m = max(mo,mi), then the condition of the theorem 
becomes mi > m — (d/2); in the case H = 1/2 this is known from [3]. Thus, beside 
extending the results of [3] to fractional-in-time noise, we also generalize the necessary 
and sufficient condition for consistency of the estimator. 

While parameter estimation for the finite-dimensional fractional OU and similar pro- 
cesses has been recently investigated by Tudor and Viens [H] for all H G (0, 1), our 
analysis in infinite dimensions requires more delicate results: an explicit expression 
for the Laplace transform of a certain functional of the fractional OU process, as 
obtained by Kleptsyna and Le Brenton [6], and for now this expression exists only 
for H > 1/2. 



2. Stochastic Parabolic Equations with Additive FBM 

In this section we introduce a diagonalizable stochastic parabolic equation depending 
on a parameter and study the main properties of the solution. 

Let H be a separable Hilbert space with the inner product (■, -)o an d the corresponding 
norm || ■ || . Let A be a densely-defined linear operator on H with the following 
property: there exists a positive number c such that ||Am|| > c||«||o for every u from 
the domain of A. Then the operator powers A 7 , 7 G R, are well defined and generate 
the spaces H 7 : for 7 > 0, H 7 is the domain of A 7 ; H° = H; for 7 < 0, H 7 is the 
completion of H with respect to the norm || ■ || 7 := ||A 7 ■ || (see for instance Krein 
at al. [8]). By construction, the collection of spaces {H 7 , 7 e R} has the following 
properties: 

o A 7 (H r ) = H r " 7 for every 7, r E R; 

o For 7! < 72 the space H 72 is densely and continuously embedded into H 71 : 

H 72 C H 71 and there exists a positive number cyi such that ||w|| 7l < Ci2||w|| 72 

for all u G H 72 ; 
o For every 7 G R and m > 0, the space H 7_m is the dual of H 7+m relative to 

the inner product in H 7 , with duality (■, -} Tjm given by 



(ui,U2) 7 ,m = (A 7 - m Ml ,A 7+m M 2 ) , where u x G H 7 ~ m , u 2 G H 



7+m 
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In the above construction, the operator A can be bounded, and then the norms in all 
the spaces H 7 will be equivalent. A more interesting situation is therefore when A is 
unbounded and plays the role of the first-order operator. 

Let (Q, J 7 , P) be a probability space and let {wf, j > 1} be a collection of independent 
fractional Brownian motions on this space with the same Hurst parameter H G (0, 1): 

Ewf(t) = 0, E(wf(t)wf(s)) = \ (t 2H + s 2H -\t- s\ 2H ) . 

Consider the following equation: 

' du(t) + (A) + 0Ai)u{t)dt = J2 9j(t)dwf(t), < t < T, 

j>i (2.1) 

^w(O) = u 

where Aq, Ax are linear operators, gj are non-random, and 9 is a scalar parameter 
belonging to an open set 0Cl. 

Definition 2.1. 

(1) Equation (12.11) is called diagonal izable if the operators Aq, Ax, have a 
common system of eigenf unctions {hj, j > 1} such that {hj, j > 1} is an 
orthonormal basis in H and each hj belongs to f] R H 7 . 

(2) Equation f l2.ll) is called (m, 7) -parabolic for some numbers m > and'-/ G M 

if 

o the operator Aq + 6A\ is uniformly bounded from H 7+m to H 7 m for 
9 G © : there exists a positive real number C\ such that 

\\(Ao + 9Ai)v\\^ m < Ci|H| 7+m (2.2) 

for all9e&,ve H 7+m ; 
o there exists a positive number S and a real number C such that, for every 
v g H 7+ "\ 9eQ, 

- 2((A + 9A 1 )v, v) 7 , m + 5\\v\\ 2 J+m < C\\v\\ 2 r (2.3) 

Remark 2.2. If equation (12.11) is (m, 7)-parabolic, then condition (12.31) implies that 

({2Ao + 29 Ax + CI)v, v) 7 , m > 6\\v\\ J+m , 

where / is the identity operator. The Cauchy-Schwartz inequality and the continuous 
embedding of H 7+m into H 7 then imply 

\\(2A + 29Ai + C/)i;|| 7 > 5i||v|| 7 

for some b~i > uniformly in 9 G O. As a result, we can take A = (2Ao + 29* Ai + 
(jjy/( 2m ) f or gome fixed 9* G 0. If the operator Ao + 9A± is unbounded, it is natural 
to say that Ao + 9A\ has order 2m and A has order 1. 

From now on, if equation (12. ip is (m,^) -parabolic and diagonalizable, we will assume 
that the operator A has the same eigenf unctions as the operators Ao, Ax; by Remark 
\2.2\ this leads to no loss of generality. 

For a diagonalizable equation, condition (12. 3p can be expressed in terms of the eigen- 
values of the operators in the equation. 
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Theorem 2.3. Assume that equation ( 12.11) is diagonalizable and 

Aohj = Pjhj, A-ihj = Vjhj. 

With no loss of generality (see Remark I2.2[) . we also assume that 

Ahj = Xjhj. 

Then equation (12.11) is (m, 7) -parabolic if and only if there exist positive real numbers 
5, C\ and a real number C<i such that, for all j > 1 and &£0, 

\j 2m \Pj + 0"j\ < Ci; (2.4) 

- 2{p 5 + 9vj) + 5X 2m < C 2 . (2.5) 

Proof. We show that, for a diagonalizable equation, (12.41) is equivalent to (12. 2p and 
(12.51) is equivalent to (12.31) . Indeed, note that for every 7, r6l, 

Then fl2l} is (Q and (J22D is (Q, with v = hj. Since both flUD and ([23]) are 
uniform in j and the collection {hj, j > 1} is dense in every H 7 , the proof of the 
theorem is complete. □ 

Remark 2.4. (a) As conditions (j2.4p . (12.51) do not involve 7, we conclude that a 
diagonalizable equation is (m, 7)-parabolic for some 7 if and only if it is (m, 7)- 
parabolic /or every 7. As a result, in the future we will simply say that the equation 
is m-parabolic. 

(b) If the operators Aq + 6A\ and A are unbounded, then ( 12.51) implies that Hj{6) = 
Pj + 0Uj is positive for all sufficiently large j. 

From now on we will assume that equation (12. ip is diagonalizable and fix the basis 
{hj, j > 1} in H. Since each hj belongs to every H 7 and, by construction, f] H 7 is 
dense in M H 7 , every element / of [LH 7 has a unique expansion J2j>i fjhj, where 
fj = (f, hj)o, m f° r a suitable m. 

Definition 2.5. The space-time fractional Brownian motion W H is an element 
of U 7 gK H 7 with the expansion 

W H (t) = ^f(t)h r (2.6) 

Definition 2.6. Let W H be a space-time fractional Brownian motion. The solution 
of the diagonalizable equation 

' du(t) + {Aq + 6Ai)u(t)dt = dW H {t), < t < T, 
u(0) = u 

uq G H, is a random process with values in M H 7 and an expansion 

u{t) = Y,^m v (2.8) 
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where 

Uj (t) = («o, h j ) e- ( - e ^ +p ^ t + / e-^+nW-^dw? (s). (2.9) 



Notice that, due to the special structure of the equation, Definition 12.61 implies both 
existence and uniqueness of the solution. 

To simplify further notations we write 

li j {0) = 0v i + Pj . (2.10) 

By (12. 5p . if equation (12.11) is m-parabolic and diagonalizable, then, for every 6 G 0, 
there exists a positive integer J such that 

fjLj(9) > for all j > J. 

Theorem 2.7. Assume that 

(1) H > 1/2; 

(2) equation (12. ip is m-parabolic and diagonalizable; 

(3) There exists a positive real number 7 such that 

£(l + |/i^)|)-^<oo. (2.11) 

Then, for every t > 0, 

(1) W H (t) eL 2 (Q;B.- m -y); 

(2) u(t) eL 2 (^H- m7+2mH ). 

Proof. Condition (12.111) implies that lim^oo \fij\ = 00, and consequently the operators 
.4o + #-<4.i and A are unbounded. The parabolicity assumption and Theorem 12.31 then 
imply that, for all sufficiently large j, 

1 + |^(0)| <C 2 Af\ 
uniformly in 9 G 0. 

E||^(t)|| 2 _ m7 = t™J2 A" 2 " 17 < ^£(1 + I^WI)" 7 < °°- 

Next, the properties of the fractional Brownian motion imply 

Eu%t)=H(2H-l)e- 2 ^ e)t f f e^ {s ^\ Sl - s 2 \ 2H - 2 d Sl ds 2] 

Jo Jo 

see, for example, Pipiras and Taqqu [TTJ, formulas (4.1), (4.2)]. By direct computation, 
lim \nj[e)\ 2Hl M.u 2 At) = H(2H - 1) / x 2H - 2 e~ x dx = H(2H - 1)T(2H - 1). (2.12) 

Consequently, 

00 

J2( l + \^mr i+2H nuAt)\ 2 <oo, (2.13) 

and the second conclusion of the theorem follows. D 
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Example 2.8. (a) For < t < T and x G (0, 1), consider the equation 

du{t,x) -9u xx (t,x)dt = dW H (t 1 x) (2.14) 

with periodic boundary conditions, where u xx = d 2 u/dx 2 . Then H 7 is the Sobolev 
space on the unit circle (see, for example, Shubin [13j Section 1.7]) and A = \/I — A, 
where A is the Laplace operator on (0, 1) with periodic boundary conditions. Direct 
computations show that equation (12.141) is diagonalizable; it is 1-parabolic if and only 
if > 0. Also, fjLj = —9n 2 j 2 , so that, by Theorem 12 .7! the solution u(t) of f)2.14|) is an 
element of L 2 (fi; H~ 7+2H ) for every t > 0, 7 > 1/2, and > 0. 

(b) Let G be a smooth bounded domain in R d . Let A be the Laplace operator on G 
with zero boundary conditions. It is known (for example, from Shubin [T3j ) , that 

(1) the eigenfunctions {hj, j > 1} of A are smooth in G and form an orthonormal 
basis in L 2 {G); 

(2) the corresponding eigenvalues <jj, j > 1, can be arranged so that < —C\ < 
— o"2 < . . ., and there exists a number c > such that \<jj\ ~ cj 2 ^ d , that is, 

lim I ctj I r 2/d = c. (2.15) 

We take H = L 2 (G), A = yl — A, where I is the identity operator. Then ||Aw|| > 
y/1 — o\\\u\\o and the operator A generates the Hilbert spaces H 7 , and, for every 
7 G R, the space H 7 is the closure of the set of smooth compactly supported function 
on G with respect to the norm 

, V2 



o>i 



1 + j 2 ) J \(pj\ 2 I , where (fj — / (p(x)hj(x)dx, 



which is an equivalent norm in H 7 . Then, for every 9 G R, the stochastic equation 

du- (Au + 9u)dt = dW H (t,x) (2.16) 

is diagonalizable and 1-parabolic. Indeed, we have At — I, Ao — —A, and 

-2(AoV,v) %1 = -2\\v\\* +1 +2\\u\\*, 

so that O holds with 6 = 2 and C = 2 - 9. Finally, by (12351) we see that (12111 
holds for every 7 > d/2. As a result, by Theorem 12.71 the solution u(t) of (12.161) is 
an element of L 2 (fi; H" 7+2H ) for every t > 0, 7 > d/2, and G R. 

3. The Maximum Likelihood Estimator and its Properties 

Consider the diagonalizable equation 

du(t) + (A + 9A l )u[t)dt = dW H [t) (3.1) 

with solution u{t) = J2j>i u j(t)hj gi ven by (12.91) ; for simplicity, we assume that 
w(0) = 0. Suppose that the processes Ui(t), . . . , Ujsr(t) can be observed for all t G [0, T]. 
The problem is to estimate the parameter using these observations. 
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Recall the notation Pj(9) = Pj + Vj9, where pj and Uj are the eigenvalues of Aq and 
.4.1, respectively. Then each Uj is a fractional Ornstein-Uhlenbeck process satisfying 

d Uj {t) = -p^u^dt + dwf{t), Uj {0) = 0, (3.2) 

and, because of the independence of wf for different j, the processes U\, . . . , Un are 
(statistically) independent. 

Let T denote the Gamma- function (see (j2.12p ). Following Kleptsyna and Le Brenton 
[6], we introduce the notations 

k h = 2HT {^-h)t(h+^Y k H {t, s) = K H \s^- H (t - s^- H ; (3.3) 

A " = r(g-V) ' w " (t) = H ; ( } 

Mf{t) = j k H (t,s)dwf(s), Q 3 (t) = ^ f k H (t,s) Uj (s)ds; (3.5) 

Zj(t) = / k H (t,s)duj(s). (3.6) 

Jo 

By a Girsanov-type formula (see, for example, Kleptsyna et al. [7J, Theorem 3]), the 
measure in the space of continuous, M^-valued functions, generated by the process 
(ui, . . . ,un) is absolutely continuous with respect to the measure generated by the 
process (w± , . . . , w§), and the density is 

e *P f-f>(0) [ QAs)dZAs)-j^^^ f\){s)d^ H {sA . (3.7) 

Maximizing this density with respect to 9 gives the Maximum Likelihood Estimator 
(MLE): 

N pT 



Yl / "iQiOO ( dZ A s ) + PjQj(s)dw H (s)) 

7 = 1 Jo 



'JV 



N 



(3.J 



i=i ■ /o 

An important feature of (13. 8p is that the process Zj is a semi-martingale (P, Lemma 
2.1]), and so there is no stochastic integration with respect to fractional Brownian 
motion: f VjQj(s)dZj(s) is an Ito integral. Notice that, when H = 1/2, we have 
kg = 1, Wjy(s) = s, Qj(s) = Zj(s) = Uj(s), and (13. 8p becomes 



TV r 

X / "j%00 (d«j(s) + p j u j (s)ds) 

_•_ i <>0 



'TV 



,=1^0 



TV T 



(3.9) 



which is the MLE from [3]. 
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Let us also emphasize that an implementation of ( 13.81) is impossible without the 
knowledge of H . 

The following is the main result of the paper. 



Theorem 3.1. Under the assumptions of Theorem 2.7, the following conditions are 
equivalent: 

(2) lim 8n = with probability one, (3.11) 

where J = min{j : Hi{6) > for all i > j}. 

Proof. Following Kleptsyna and Le Brenton [SJ Equation (4.1)], we conclude that 

h - » = - £ H^ WM (3-12) 

Both the top and the bottom on the right-hand side of (13.121) are sums of independent 
random variables; moreover, it is known from [6] page 242] that 

Elf Q 3 {s)dMf{s)\ =E f Q 2 {s)dw H {s)ds. (3.13) 

From the expression for the Laplace transform of J Q Q 2 (s)dw H (s)ds (see [6] Equation 
(4.2)]) direct computations show that 

lim //,-(0)E / Q 2 Js)dw H {s)ds = ^ > (3.14) 

j^oo J 2 

and, with Var(£) denoting the variance of the random variable £, 

lim /4(fl)Var ( / Q 2 (s)dw H (s)ds] = - > 0; (3.15) 

i^oo \Jo / 2 

a detailed derivation of (13.141) and (I3.15P is given in the appendix, Lemmas IA.1I and 
IA.2I respectively. 

We now see that if (13.101) does not hold, then, by (I3.14p . the series 

J2 f ^Q%s)dw H {s) 

3>1 J ° 

converges with probability one, which, by ( 13.121) . means that ( 13.111) cannot hold. 
On the other hand, if (13.101) holds, then 

— ,2,,-l N 
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Indeed, setting a n = ^/i^ 1 and A n = YTj=i a j> we notice that 

Al ~ ^ \A n An-! J ~ A,' 
n>J n n>J+l V n n ly J 

Then the strong law of large numbers, together with the observation 

E / Q J (s)dMf(s) = 0, j>l, 
Jo 

implies 



Eli Jo VjQ^dMjis) 

lim — -rf ~ = with probability one. 



Next, it follows from (13TT6]) and (12TTD that 

,4, ,-3 



n>J £ ^,V 



-j=J 

Then another application of the strong law of large numbers implies that 

lim E^J^M^) = t (3 , 8) 

with probability one, and (13. lip follows. D 

Corollary 3.2. Under assumptions of Theorem \2.1\ if (13.101) holds, then 



lim 

AT— »oo 



N o 

vi 



m distribution, where ( is a Gaussian random variable with mean zero. 

Proof. This follows from ( 13.121) . (13.181) . and the central limit theorem for the sum of 
independent random variables. □ 

Let us now consider a more general equation 

du = (A + 9Ai)udt + BdW H (t), 

where B is a linear operator. If B~ x exists, the equation reduced to (13. ip by considering 
v = B^ 1 u. If B~ x does not exist, we have two possibilities: 

(1) (u , hi)o = for every i such that Bh{ = 0. In this case, Ui(i) = for all t > 0, 
so that we can factor out the kernel of B and reduce the problem to invertible 
B. 

(2) («o, hi)o 7^ for some i such that Bhi = 0. In this case, Ui(t) = Ui(0)e~ Pit ~ Uldt 
and 9 is determined exactly from the observations of Ui(t): 

a 1 i M *( g ) Pi ,, 
6 = — -, r In — -r-r- , if s. 

Viit-S) Ui{t) Vi 
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Let Aq, A\ be differential or pseudo-differential operators, either on a smooth bounded 
domain in M. d or on a smooth compact d- dimensional manifold, and let mo, mi, be 
the orders of Aq, A\ respectively, so that 2m = ma.x(m ,mi). Then, under rather 
general conditions we have 

lim \vA] mi/d = Cl , lim ^{9)j 2m/d = c{9) (3.20) 

for some positive numbers c\,c{9); see, for example, Il'in [5] or Safarov and Vassiliev 
[12J. In particular, this is the case for the operators in equations (12.141) and (j2.16p . 

If (I3.20p holds, then condition (13.101) becomes 

m l >m-(d/2), (3.21) 

which, in the case H = 1/2, was established by Huebner and Rozovskii [3j. In 
particular, (I3.2ip holds for equation (12.141) (where 2m = m\ = 2), and for equation 
fl2~T6|) if d > 2 (where 2m = 2,m x = 0). 

Note that, at least as long as H > 1/2, conditions (13.101) and (13.211) do not involve 
H. 

The maximum likelihood estimator (13. 8p has three features that are clearly attractive: 
consistency, asymptotic normality, and absence of stochastic integration with respect 
to fractional Brownian motion. On the other hand, actual implementation of 



is problematic: when H > 1/2, computing the processes Qj and Zj is certainly 
nontrivial. Estimator (13.91) is defined for all H > 1/2 and contains only the processes 
Uj, but, when H > 1/2, is not an MLE and is even harder to implement because of 
the stochastic integral with respect to Uj. 

With or without condition (13.101) . a consistent estimator of 6 is possible in the large 
time asymptotic: for every j > 1, 

Um Jo VjQj( s ) i dZ A s ) + PjQj(s)dw H (s)) = _ 
T ^ j^ 2 3 Q){s)d^ H {s) 

with probability one ([61 Proposition 2.2]). For H > 1/2, implementation of this 
estimator is essentially equivalent to the implementation of (13. 8p . 

An alternative to ( 13.22J) was suggested by Maslowski and Pospisil [10] using the er- 
godic properties of the OU process. Let us first illustrate the idea on a simple example. 

If a > and w — w(t) is a standard one-dimensional Brownian motion, then the OU 
process dX = —aX(t)dt + dw(t) is ergodic and its unique invariant distribution is 
normal with zero mean and variance (2a) _1 . In particular, 



r 
with probability one, and so 



lim i f X 2 {t)dt = ±- (3.23) 



a(T) = — ^- (3.24) 

2j T X%t)dt 
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is a consistent estimator of a in the long-time asymptotic. Note that the maximum 
likelihood estimator in this case is 

, s [J X{t)dX(t) 

= _ Jo w w 3 25 

S*X\s)ds 

and is strongly consistent for every o 6 1 [9, Theorem 17.4]. 

Similarly, if a > 0, then the fractional OU process 

dX(t) = -aX{t)dt + dw H (t), X(0) = (3.26) 

is Gaussian, and, by (I2.12p on page El converges in distribution, as t — > oo, to the 
Gaussian random variable with zero mean and variance c(H)a~ 2H , where 

c(H) = H(2H-l)T(2H-l); (3.27) 

notice that, in the limit H \ 1/2, we recover the result for the usual OU process. 
Further investigation shows that, similar to (I3.23p . 

(see [10j). As a result, for every j such that 0Uj + pj > 0, we have 

1 f T 9/ n , c(H) 

with probability one. Under an additional assumption that Vj ^ 0, we get an estima- 
tor of 9 

0U)( T ) = 1 ( c(iJ)r ^1 2H - Pi (3 29) 

This estimator is strongly consistent in the long time asymptotic: Ymit^oo \0^'{T) — 
^|=0 with probability one ( [TU1 Theorem 5.2]). While not a maximum likelihood 
estimator, (I3.29P is easier to implement computationally than f)3.8p . If, in Theorem 
12.71 on page El we have A® = 0, Uj > 0, and 7 < 2H, then a version of f)3.30p exists 
using all the Fourier coefficients Uj, j > 1 : 

see [TUl Theorem 5.2]. 

An interesting open question related to both f)3.8p and f)3.29p . (13.301) is how to combine 
estimation of 9 with estimation of H. 
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Appendix 
Below, we prove equalities (13.141) and (13.151) . 
Lemma A.l. For every 9 G O and H G [1/2, 1), 

lim ^.(0)1 / Q){s)dw H {s) 



j^oo Jo 2 

r-T 



Proof. Denote by \l/y(a,/ij) the Laplace transform of J Q Qj(s)dw H (s), namely 

V$(a,n J (9)) = 'Eexpl-a Q 2 j (s)dw H (s) i , a > 0. (A.l) 

We will use the expression for \l/^ from [61 page 242], and write it as follows 



(^,-a)T 



\l/y (a, /j,j) = ae 2 [A^(/ij,a 
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where pij = fij(9), a := *l jjq + 2a 



Af(^,a) 



iraTe (a — \x 



aT 



■■I-„\^)I„-A^ 



aT 



+ e 



-aT 



4sm(irH) 

. , (aT\ , /aT\ 

a sinh I — I + Hj cosh I — I 



and I p is the modified Bessel function of the first kind and order p. 

Note that 



E / Q;;( s )dwtf( s ) 



d^(a,/j,j) 



da 



a=0 



Direct evaluations (for example, using Mathematica computer algebra system) give 



da 



a=0 



4/x?e^ 



where csc(x) = l/sin(x). By combining formulas (6.106), (6.155), and (6.162) in pQ, 
we conclude that, for all p E (— 1, 1), p ^ 0, we have I p (x) ~ e x /\j2nx, x — > oo, that 
is, 



lim V27TX e x /«(a;) = 1. 

x— >+oo 



(A.2) 



Therefore 










d^§(a,fij) 


2 + 2e^ T (l - ///T) - 


e^ T 


csc(if7r) 


T 


da 


a=o 4/i|e« T 






2/ij 






a=0 


T 

= ~2' 





J ^oo, 



and the lemma is proved. 

Lemma A.2. For every 9 E and H E [1/2, 1) 



a 



J^OO 



T 



lim /iJ(0)Var / Q?(s)dw ff (s; 



T 

2" 



Proof. Note that 

V:=Vax([ Q^dwffis)) 



da 2 



<h, 



d 2 ^{a,Hj) (dV$(a,(ij) 



da 



(A.3) 



a=0 
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with \1/^ from ( 1A.1I) . Direct evaluation of the right hand side of ( 1A.3I) (for example, 
using Mathematica computer algebra system) gives 

V = pi^- (2 - 8e"' T (l + N T) + 2e 2 ^(-5 + 2//,T) 

+ m TcBc{«H) [ - 2e^ T ^Th_ H (^f\ I E _ X (^f\ 
+ I-e i^f) {4(-l + e" r (l + n/r))I B -i (^) 

" 2^2^ (^f ) + ^T/^ (^) I_ fl (^) csc^tt)}]), 

where csc(a;) = l/sin(a;) and I p is the modified Bessel function of the first kind and 
order p. 

Using ( 1A.2I) . we conclude that 



3 _ 3 /-10 + 4csc(if7r) + csc 2 (Hn) 



lim fi 6 j (9)Y = lim ^ ( — ■ K -—{ — + 

+ 



csc{Hir) + 2 + 2/ijT T 



2//fe^ T ' 2/x 3 



_T 
~"2 
and complete the proof of the lemma. □ 
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